Scaling of entanglement entropy across Lifshitz transitions 
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We investigate the scaling of the bipartite entanglement entropy across Lifshitz quantum phase 
transitions, where the topology of the Fermi surface changes without any changes in symmetry. We 
present both numerical and analytical results which show that Lifshitz transitions are characterized 
by a well-defined set of critical exponents for the entanglement entropy near the phase transition. 
In one dimension, we show that the entanglement entropy exhibits a length scale that diverges as 
the system approaches a Lifshitz critical point. In two dimensions, the leading and sub-leading 
coefficients of the scaling of entanglement entropy show distinct power-law singularities at critical 
points. The effect of weak interactions is investigated using the density matrix renormalization 
group algorithm. 
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I. INTRODUCTION 

Quantum phase transitions are characterized by quali- 
tative changes in the nature of the ground state as a pa- 
rameter in the Hamiltonian is tuned across the transition, 
which mathematically manifests as non-analyticities in 
the ground-state energy and other local operators in the 
thermodynamic limit [1, 2]. The conventional framework 
for studying different phases of matter focuses on dif- 
ferent symmetries and the identification of local order 
parameters that spontaneously break them. However, 
the discovery of distinct phases of matter with the same 
symmetry has recently introduced a more general notion 
of order associated with the global topology of ground- 
state wave functions [3]. For example, different quantum 
Hall states [4, 5] and topological insulators [6] are charac- 
terized by topological invariants. Topologically ordered 
states have distinct patterns of long range entanglement 
which can be probed using various measures of entangle- 
ment. In particular, one can consider a subsystem of 
the full system with ground-state wave function | <L) and 
its associated reduced density matrix pn = Tr^d*) 1 ) 
obtained by tracing out the degrees of freedom in the re- 
mainder of the system, f2. The bipartite entanglement 
entropy is defined as the von Neumann entropy of the 
reduced density matrix, 

5(0) = -Tr^ln fa), (1) 

and is symmetric (for the pure state |\E')) between Q and 
ft. Non-trivial topological order can be then detected 
from the universal component of the entanglement en- 
tropy [7-10]. 

As entanglement entropy can serve as a useful di- 
agnostic for non-trivial gapped phases, it can also be 
used to detect non-trivial topological structures in gap- 
less phases. Many interesting properties of the entangle- 
ment entropy, including striking universal scaling behav- 
ior in one-dimensional conformal field theories (CFTs), 



have been investigated over the years [11, 12]. More 
recently, entanglement entropy has been computed for 
various gaplcss states in higher dimensions [13-18]. In 
this work, we use entanglement entropy to detect and 
characterize quantum phase transitions between two gap- 
less phases with the same symmetry but with different 
topologies. In particular, we study the scaling of the bi- 
partite entanglement entropy across Lifshitz transitions, 
an example of a quantum phase transition in which there 
is no change in symmetry [19-21]. Instead, Lifshitz tran- 
sitions arc driven by changes in the topology of the Fermi 
surface, for example when two Fermi surfaces merge to 
form a single Fermi surface. The transition is signaled by 
van Hove singularities in the density of states and is due 
to the critical points where the Fermi velocity de(k)/dk 
vanishes at the Fermi surface. We give both analytical 
and numerical results for the scaling of entanglement en- 
tropy in some simple systems that possess Lifshitz tran- 
sitions and show that the entanglement entropy exhibits 
critical behavior across such transitions. We also use 
the density matrix renormalization group (DMRG) algo- 
rithm to investigate the effect of weak interactions. 

It is worth noting that, because of the equivalence be- 
tween entanglement entropy and particle number fluc- 
tuations in free-fermion systems [22-24], many of the 
conclusions regarding entanglement entropy in the sys- 
tems studied below apply to particle number fluctua- 
tions. In addition, recent proposals for measuring the 
entanglement entropy in interacting systems [25-28] of- 
fer a promising way to observe the effects investigated in 
this work. 



II. LIFSHITZ TRANSITIONS IN ID 

We begin with the simple case of one dimension. Con- 
sider the following, infinite system of one-dimensional 
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FIG. 1: (color online) Occupied modes (shaded) for the dispersion corresponding to Eq. (2) with n = 2. 
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FIG. 2: (color online) Occupied modes (shaded) for the dispersion corresponding to Eq. (2) with n = 3. 



fcrmions 



[(cJci+l+*4ci+n)+ h -C-], 



(2) 



where cj , Cj are fermionic creation and annihilation op- 
erators on site i and t is the amplitude of the n-th 
nearest-neighbor tunneling relative to the amplitude of 
nearest-neighbor tunneling. The system has the disper- 
sion relation e n (k) = — 2[cosfc + tcos(nk)]. Since the 
dispersion is symmetric about k = 0, let the occupied 
modes in the ground state be k lying in the positive in- 
tervals (0, k\), (k2, £3), . . . , (k m -i, k m ) if m is odd and 
(0, fci), (k 2 , fc 3 ), . . . (fc m _ 2 , fcm-i), (k m ,ir) if m is even, 
and similarly for negative values of k. Here we assume 
t > 0. As shown for the case of n = 2 in Fig. 1 and 
n = 3 in Fig. 2, m is the number of Fermi surface "com- 
ponents" in the Brillouin zone, i.e., the Fermi surface has 
2m points. 

The entanglement entropy in the ground state of a sys- 



tem of non-interacting fcrmions can be determined from 
the matrix of Green's functions [29] 



MijiSl) = (etc,-) - 



2tt 



dk Q(-e n {k))e 



-ik(i-j) 



(3) 



for i, j e n, where f2 is the subsystem of interest and 
consists of an interval of length £ (lattice sites). Here 
9(a;) is the Heaviside step function. The entanglement 
entropy is then given by 

S(£) = — Tr[MlnM + (1 - M)ln(l - M)\. (4) 

It is a remarkable fact that Eq. (4) can be evaluated an- 
alytically in the limit of large £, using techniques from 
the theory of Tocplitz, or translationally-invariant, ma- 
trices, in particular the Fisher-Hartwig conjecture [30]. A 
straightforward generalization of the calculation for the 
case t = carried out in Ref. [30] reveals that for t > 
and i!^oowe have 



S(£) 
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where s\ = (ln2)/3 + Ti with 



Ti 
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dt 



1 cosh(t/2) 
3i tsinh 2 (t/2) 2sinh 3 (t/2) 



+ 



• (6) 



The 
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-l , 



related to the Fermi momenta by 
cos~" Xj where \xj\ < 1 are the roots of 
tT n (x) + x = with T n (x) the Chcbyshev polynomials 
of the first kind. Eq. (5) is valid if £ is sufficiently large 
relative to \t — t c \, the distance from the critical point. 
The main feature is the discontinuity of the logarithmic 
term: the prcfactor of the leading, logarithmic term is 
proportional to the number of Fermi surface components 
to. That this is the special case of a more general formula 
will become clear when we consider the Gioev-Klich 
formula in the next section. 



which is consistent with Eq. (5) where the prefactor of the 
logarithmic scaling was derived through a microscopic 
calculation. 

In two dimensions Eq. (8) can also be simplified. For 
convenience we consider an I x £ square region of real 
space so that the scaled region Q is a unit square. The 
normal unit vectors in real space are thus ±x, ±y. Since 
the Fermi surface is a curve in two dimensions, we can 
parametrize the curve by k x = k x {0) and k y — k y (6), 
with the unit normal vector given by the normalized ve- 
locity rife = Vk£/|Vk£| where e(k) is the dispersion rela- 
tion. Then 



Jan 



1 + 



dA x |n fc ■ n x | = 
Meanwhile, the line element is given by 



(10) 



III. THE GIOEV-KLICH FORMULA IN ID AND 

2D 

Various arguments [31-33] have demonstrated that 
the leading scaling of the entanglement entropy in free 
fcrmions with subsystem size exhibits a logarithmic vio- 
lation of the area law with subsystem size (linear dimen- 
sion) £, S(£) <~ £ d ~ 1 \n£ in d dimensions. More specif- 
ically, a conjecture by Widom [34] leads to the explicit 
formula [32] 

S(£) ~ a£ d - 1 \n£, (7) 
a = ( old-i To I dAk I dA x \n k -h x \, (8) 

for the leading behavior of the entanglement entropy in 
an infinite system of gapless fermions with a Fermi sur- 
face of co-dimension 1, i.e., the surface separating the 
occupied and unoccupied regions of the rf-dimensional 
Brillouoin zone is a {d— l)-dimensional surface. Here dT 
is the Fermi surface, <9fi is the real-space boundary of 
the subsystem fi scaled such that the volume of £1 is 1, 
and fife and h x are unit vectors normal to dT and OQ, 
respectively. Although the mathematical derivation of 
Eq. (8) is somewhat complex, an intuitive explanation 
was given in Ref. [33] by associating each patch of the 
(d — l)-dimensional Fermi surface with a chiral confor- 
mal field theory that contributes (1/6) In £ to the total 
entanglement entropy, so that a simple patch-counting 
argument leads to the result. 

In one dimension Eq. (8) is simple to evaluate. For 
a real-space interval of length £ embedded in an infinite 
system, the geometric factor J gr dA k J dQ dA x \h k ■ ii x \ 
counts the number of boundary points on the Fermi sur- 
face and the real-space boundary, or 4m where m is the 
number of Fermi surface components. Therefore 



m 

am = j, 



(9) 



but since 



dS k = 



de 
00 



' dk x 

de 



de dk x 
dk x dO 



dk y 
d9 



de dk y 
dky d6 



<W, 



along the Fermi surface, we have 

de/dky dk x /de 



de/dk x 



and we can write 
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dfx r 



de 



de 



(11) 



(12) 



(13) 



(14) 



Eq. (14) emphasizes the purely geometrical character of 
the prefactor of the leading scaling of the entanglement 
entropy in such systems, and can be evaluated by split- 
ting the parametrization of the Fermi surface into piece- 
wise monotonic regions. As the simplest example, con- 
sider the case of nearest-neighbor hopping on the square 
lattice at zero chemical potential. Then the dispersion 
is e(k) = — 2{cosk x +cosk y ). Since the Fermi surface 
is invariant under 90-degree rotations (which we will as- 
sume in all of the models considered in this work), we can 
treat the first quadrant only and multiply by 4. In the 
first quadrant of the Brillouin zone the Fermi surface is 



clearly parametrized by k x 
so that Eq. (14) gives 



k y = 7r - e for 6 e [0,7r), 
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k x {e) 



k v {6) 
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3' 



(15) 



More generally, including a chemical potential term in 
the dispersion leads to the prefactor [35, 36] 
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7T 



(16) 
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FIG. 3: (color online) Fermi surfaces (blue, thick lines) for the dispersion in Eq. (17) with fj, = 0. Solid (dashed) lines represent 
positive (negative) energy states. 



for < 4, which has a cusp at fi = such that near 
/i = 0we have a(fi) w a(0) - (2/37r)|/i| 1 / 2 . Thus at 
fi = we have a simple example of a Lifshitz transition, 
and in fact this transition is very similar in character to 
the transition at t — studied below, cf. Figs. 3(c), 3(d). 
Consider now the more interesting dispersion relation 



e(k) = — 2(cosfc a; +cosfc y ) 

— 2t[cos(2k x ) + cos(2fc y ) 



fi, 



(17) 



which arises from nearest-neighbor and next-nexi-nearest 
neighbor coupling on the square lattice. We treat this 
case to simplify the algebra, although a slightly more 
realistic model would include next-nearest neighbor cou- 
pling instead. To parametrize the Fermi surface, we again 
use the 90-degree rotational symmetry to limit ourselves 
to the first quadrant of the Brillouin zone and multiply 
the overall result by 4. Let (we assume f ^ for the 
moment) 



w u 



R{t, fx) = y/i + not) + 2C(ty 

Then the Fermi surface is defined by 



k x (9) 
k y {9) 



cos _1 [C(t) +R(t,fi) cosfl], 
cos _1 [C(i)+ii(i,/z) sin 6] 



(18) 
(19) 



(20) 
(21) 



with the ranges of 9 depending on t,fi so that —1 < 
C + R cos 9 < 1 and — 1 < C + R sm6 < 1. For example, 
if fi = then 



G < 



' (vr/2 + W(t), 2tt - W(t)) and 

(W(t), %/2-W(t)) 
(n/2 + W(t), 2ir-W(t)) 
(-n/2 + W(t), n-W(t)) 
(-tt/2 + W(t), 7T - W{i)) and 
(-7r + W(t), -7r/2-W(t)) 



for t < -1, 
for -1 < t < 0, 
for < f < 1, 



for t > 1, 



where 



= cos" 



1 



+ C(I*I) 
R(\t\) 



(22) 



(23) 



We note that the ranges for t > are merely shifted by 
7r from the values for t < 0. Fig. 3 shows some represen- 
tative Fermi surfaces for different values of t. For t < — 1 
the occupied modes contain a "hole" that disappears at 
t = — 1, while at t = 1 disconnected Fermi components 
appear. 

The result of using Eq. (14) is invariant under t — > — t 
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FIG. 4: The coefficient of the leading contribution to the 
entanglement entropy, given by Eq. (24). Dashed lines show 
the square-root nature of the kinks/cusps. 



and can be described as 

(l-f(W(t),t) for < \t\ < 1/4, 

1 - f(W(t), t) + 2f(W(t),t) - 2/(tt/2, t) 
for 1/4 < \t\ < 1, 
1 - f(W(t), t) + 2f(W(t),t) - 2/(7r/2, t) 
{-f(W(t)-7r,t) + l for jtj > 1, 

(24) 



a[t) = - x< 



where ki is the «-th component of momentum measured 
from the momentum point where the Lifshitz transition 
occurs at 5t = and p is an integer that parameter- 
izes the topology of the Fermi surfaces involved in the 
transition. Locally in momentum space, the dispersion 
respects SO(p, d — p) symmetry. For p = (p = d), this 
dispersion describes the transition where a new hole-like 
(electron-like) Fermi surface of topology S l<i_1 appears as 
5t becomes positive (negative). For < p < d, it de- 
scribes the transition where the neck of the hyperbolic 
Fermi surface pinches to a point to change its topology. 
Near the critical point, the area of the Fermi surface in 
momentum space changes non-analytically as 



A~ |<5i| (d - 1)/2 e(±<5f) 



for p = or d, and 



A ~ \St\ 



(d-l)/2 



(27) 



(28) 



where 



for < p < d. In d = 1, it becomes a jump as we saw 
in Sec. II. It is of note that the non-analyticity does not 
depend on p. Since the entanglement entropy is related 
to the area of the Fermi surface (specifically, each patch 
of Fermi surface contributes a logarithm [33]), the same 
non-analyticity shows up in the scaling of entanglement 
entropy. 



V. NUMERICAL RESULTS 



A. ID 



f(0,t) = - cos- 1 [C (\t\) + R(\t\) sin9}. (25) 

7T 

As shown in Fig. 4 the prefactor of the logarithmic scaling 
is continuous but non-monotonic even for the same sign 
of t (partly because particle number is not fixed as t 
changes), and moreover has square-root singularities at 
t = and \t\ = 1. Specifically, near t = we have a(t) « 
a(0) - (4/37r)|i| 1 / 2 while near \t\ > 1 we have a(t) w 
a(l) + (4/3v / 37r)(|i| - l) 1 / 2 . The presence of square-root 
kinks is a general feature in 2D, and is explained in the 
following section for arbitrary dimensions. 



IV. GENERAL DIMENSIONS 

We consider here the scaling of entanglement entropy 
near a Lifshitz transition in dimension d, focussing on 
the case where the co-dimension of the Fermi surface is 1 
and the dynamical critical exponent at the critical point 
is 2, but it is possible and straightforward to extend the 
result further to more general cases. Consider the local 
fermion dispersion 



£ (*o = 5> 2 - E 



St, 



(26) 



In this section we present results for the entanglement 
entropy computed numerically for the one-dimensional 
Hamiltonian in Eq. (2). Numerical data is obtained for 
system sizes ranging from £ ~ 20 to £ = 200 for n = 2,3 
respectively. In Fig. 5(a), we show the entanglement en- 
tropy for n = 2 as a function of t for various values of 
I. Even for relatively small £, the Lifshitz transition at 
t = 1 where the number of Fermi points changes is clearly 
visible. As £ increases, the discontinuity becomes more 
pronounced and sharper. The entanglement entropy as a 
function of £ is fit to the form a(t) \n£ + b(t), allowing us 
to extract the coefficients as a function of t. The coeffi- 
cients are shown in Fig. 5(b) and Fig. 5(c). As expected 
from the increase of the number of Fermi points by a fac- 
tor of 2, the leading coefficient a(t) doubles across the 
transition. In Figs. 6(a), 6(b), and 6(c), we show the en- 
tanglement entropy and the leading coefficients in the n 
= 3 case. Because the number of Fermi points increase 
by factor of 3, a(t) triples in magnitude across the tran- 
sition. 

Now we discuss the finite size crossover behavior of en- 
tanglement entropy near the critical point focusing on the 
case with n = 3. When t is slightly larger than the criti- 
cal hopping t c — 1/3, two small Fermi pockets are formed 
near k = ±ir/2 as is shown in Fig. 2. This sets a charac- 
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teristic length scale £* 



(t — t c ) - 1 / 2 , where fci? is 
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FIG. 5: (Color online) (a) Numerically calculated entangle- 
ment entropy in one dimension as a function of the next- 
nearest neighbor hopping t (n = 2) for various values of sub- 
system size, (b) The leading coefficient, a(t) as a function of 
t. (c) The sub-leading coefficient, b(t) as a function of t. The 
solid lines are guides to the eye. 

the size of the newly formed pockets. This can be under- 
stood from the dispersion = a(t c —t)Sk+b <5fc 3 + ..., with 
5k = k — 7r/2 near the critical point, which implies that 
the size of small Fermi pockets scale as kp ~ [t — tc) 1 / 2 . 
For I > I* , the entanglement entropy is expected to show 
the asymptotic behavior | In £, where c — 3 is the central 
charge for the low energy modes. For I < £* , however, 
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FIG. 6: (Color online) (a) Numerically calculated entangle- 
ment entropy in one dimension as a function of the next-next 
nearest neighbor hopping t (n = 3) for various values of sub- 
system size, (b) The leading coefficient, a(t), as a function of 
t. (c) The sub-leading coefficient, b(t), as a function oft. The 
solid lines are guides to the eye. 

the newly formed small Fermi surfaces are not completely 
resolved, and the entanglement entropy shows a different 
behavior. In Fig. 7, we plot the entanglement entropy 
with the constant piece subtracted, as a function of ln£ 
for various values of t with n = 3. When t is sufficiently 
close to t c , there is a wide range of I < £* where the small 
Fermi pockets are invisible so that the entanglement en- 




FIG. 7: (color online) The entanglement entropy with the 
subraction of the constant piece plotted as a function of In i 
near the critical point t c — 1/3 with n = 3. The number 
above each curve denotes the value of t. The entanglement 
entropy is shifted vertically by a different amount for each t 
to differentiate different curves. Here i ranges from 5 to 200. 
For large £, the entanglement entropy fits well with ln^ (blue 
line), whereas for small £, the entanglement entropy behaves 
as | \n£ (black line). 



tropy increases as §hi£ before it crossovers to \n£. For 
t > 0.336, the |ln^ behavior essentially disappears. In- 
terestingly, the entanglement entropy appears to increase 
with a higher power of the logarithm (\n£) A with A > 1 
at short distances before it saturates to ln^. In the non- 
interacting case with n = 3, the ground state right at the 
critical point is rather boring: it is the same as the one 
for t < t c . That is why the entanglement entropy at the 
critical point shows the usual logarithmic behavior with 
c = g. However, this is no longer the case for interact- 
ing systems. It would be of great interest in the future 
to study the possibility of the violation of the logarith- 
mic behavior in interacting quantum field theories with 
general dynamical critical exponent z > 1. 





(c) 



B. 2D 

We now turn to the two-dimensional case. The entan- 
glement entropy is calculated for the tight binding model 
on the square lattice. With the nearest-neighbor hop- 
ping whose magnitude is normalized to 1 and the next- 
next-nearest neighbor hopping t, the dispersion is given 
by Eq. (17). The chemical potential is fixed to fi = 0, 
and t is tuned to induce Lifshitz transitions. Under the 
transformation — > Cj(— the energy dispersion 
transforms as £t(fc) — » — e_t(fe). This implies that the en- 
tanglement entropy is an even function oft. We therefore 
focus on the regime with t > 0. Results are presented for 



FIG. 8: (color online) (a) Numerically calculated entangle- 
ment as a function of the next-next nearest neighbor hopping 
t for the dispersion in Eq. (17) with /i = 0. (b) The leading 
coefficient, a(t), as a function oft. (c) The sub-leading coeffi- 
cient, b(t), as a function oft. The solid line in (b) represents 
the analytical prediction from Eq. (24). 

subsystems of size £ = 5 up to £ = 80 with < t < 1.5. 
The entanglement entropy has been fitted to the general 
form 

S{£) = a(t)£]nl + b(t)£ + c{t) hit + d(t). (29) 

In Fig. 8(a), we show the entanglement entropy as a func- 
tion of t for I — 80. In Figs. 8(b) and 8(c), the first two 
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FIG. 9: (color online) Log-log plot of a(0) — a(t) (lower, blue 
line) and b(t) — b(0) (higher, red line) against t. The parameter 
t goes up to ~ 0.3. The small t behavior clearly shows a power 
law relationship for both as a function of t. The slope for the 
a(0) — a(t) line is 0.490, whereas for the bit) — 6(0) line it is 
0.335. 

leading coefficients are shown. The leading-order coef- 
ficient a(t) agree quite well with Eq. (14) as predicted 
by Widom's conjecture. In particular, a(t) — a(0) is 
found to behave as — s\t\ r near the critical point at t = 
with s « 0.40 ± 0.04 and r « 0.490 ± 0.009 as shown 
in Fig. 9, fairly consistent with the analytic predictions 
s = (4/37I-) w 0.4244 and r = 1/2. We can also extract 
the scaling behavior of the next leading-order coefficient. 
It was found that b(t) — 6(0) also follows a power law 
for t small as is shown in Fig. 9. In this case, we have 
b(t) - 6(0) ~ \t\ r with r w 0.335 ± 0.006 near t = 0. Al- 
though the sub-leading term itself is non-universal, the 
critical exponent that governs the scaling behavior near 
the Lifshitz transition may be universal. It is of great 
interest to establish this rigorously in the future. It is 
remarkable that the sub-leading coefficient vanishes with 
a smaller power than the leading coefficient. The Lifhistz 
transition at t = 1 shows similar behaviors except that 
the non-analyticity appears only for t > 1 as expected. 

We also computed the entanglement entropy for the 
model with next nearest neighbor hopping with the dis- 
persion 

e(k) = — 2(cos fc z +cos k y ) — 2t[cos(k x +k y )+coa(k x — k y )], 

corresponding to nearest-neighbor hopping. The entan- 
glement entropy and the leading coefficients are displayed 
in Figs. 10(a), 10(b), and 10(c). The behavior of t near 
the critical values is similar to that which was shown for 
the previous dispersion. 

VI. EFFECT OF INTERACTIONS 

How would the above conclusions for free fermions 
change in the presence of interactions? A priori, inter- 
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FIG. 10: (color online) (a) Numerically calculated entangle- 
ment as a function of the next nearest neighbor hopping t for 
the dispersion in Eq. (30) with /i = 0. (b) The leading coeffi- 
cient, a(t), as a function of t. (c) The sub-leading coefficient, 
b(t), as a function of t. 



actions introduce new phases that are not amenable to 
the above analysis [37-39]. In principle, the generic Lif- 
shitz critical point may disappear due to a pre-emptive 
spontaneous symmetry breaking [40]. In this case, the 
Lifshitz critical point becomes a multi-critical point that 
requires tuning more than one parameter. At sufficiently 
low energies, however, the low energy modes are still de- 
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scribed by a collection of one dimensional chiral fermions. 
In light of the argument in Ref. [33], the main features 
of the scaling of the entanglement entropy away from 
a Lifshitz transition in interacting systems with Fermi 
surfaces should be qualitatively similar for sufficiently 
large I. Indeed, the study of Lifshitz transitions in sev- 
eral mean- field approaches [21] suggests that the methods 
used in this work may be applicable to the quantitative 
analysis of such systems. It would be of interest to exam- 
ine the entanglement entropy right at the critical point 
described by interacting quantum field theories with a 
general Lifshitz scaling, z =/= 1. 




(a) V = 0.0 




. V=0.2, t=0.45 . 
0^ IT6 
l/L 



(b) V = 0.2 



FIG. 11: (color online) DMRG results for Eq. (31) with n = 
2. Results are shown for the bipartite entanglement entropy 
S(1,L) versus l/L for systems with L = 200 sites at half-filling 
with open boundary conditions, (a) Results for V = 0, (b) 
results for V = 0.2. 

In order to partly address some of these issues in one- 
dimensional case we have performed DMRG calculations 



on the model Eq. (2) in the presence of an interaction 
term V: 



H„ 



[< 



cjc. 



tic, 



) + h.c. + Vhihi 



+i 



(31) 

with hi — cjc^. A parallelized finite system DMRG 
method was used with the number of states kept ranging 
from m — 256 — 512 with open boundary conditions. By 
partitioning the system in left and right parts of linear 
extent I and L — I, respectively, the bipartite entangle- 
ment entropy S(l, L) can then be straight forwardly cal- 
culated. We first discuss our results for V = 0, n = 2 
shown in Fig. 11(a). In contrast to results in the previ- 
ous sections the DMRG calculations were performed at 
fixed particle number, at half-filling, as opposed to fixed 
chemical potential. In this case the Lifshitz transition for 
n = 2 occurs at t = 1/2. We first focus on the results for 
t = 0.45,0.49,0.50 (V = 0), which are nearly identical. 
For gapless systems with open boundaries the uniform 
part of S(l, L) is known from CFT [41] to be 



S u (l,L) = - In 
6 



"2L 




— sin 




7T 





const, 



(32) 



with c the central charge. The uniform part of our re- 
sults for t = 0.45, 0.49, 0.50 (V = 0) follow this form very 
closely with a central charge of c = 1. A known [42] 
alternating part of the entanglement entropy is also vis- 
ible in the results for t = 0.45, 0.49, 0.50 (V = 0). When 
t is increased beyond the Lifshitz transition at t = 1/2 
two effects are clearly visible: The entanglement entropy 
increases while pronounced oscillations develop. Such 
enhancement of entanglement entropy has been also ob- 
served in the spin Bose metal phase in ladder systems [43] . 
Since the system remains gapless the uniform part of 
S u (l,L) should still be given by Eq. (32), but due to the 
pronounced oscillations in S it is not possible to fit to this 
form until rather large values of t. At t = 2 such a fit 
yields an approximate central charge of c ~ 2 correspond- 
ing to a doubling of the number of Fermi points. The 
wavelength of the oscillations correspond to the momen- 
tum associated with the small Fermi surface in Fig. 1(c). 
We expect similar ocsillations to be observable in many 
other observables. Similar oscillations in the entangle- 
ment entropy have also been observed in spin chain sys- 
tems. [24, 44, 45]. 

Our results for V = 0.2 n = 2 are shown in Fig. 11(b). 
As was the case in the absence of interactions the uniform 
part of the entanglement entropy is again described by 
Eq. (32). This continues to be the case until t ~ 0.54, at 
which point the entanglement entropy again dramatically 
increases and oscillations appear. From these results for 
S(l, L) no sign of another phase appearing is apparent 
and the only effect of the interactions is to slightly shift 
the transition point away from t = 1/2. 

We have also obtained preliminary results for n = 3. In 
this case the expected tripling of the entanglement slows 
the DMRG calculations since the number of states that 
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must be kept is significantly larger. However, our pre- 
liminary results indicate similar oscillations as for n = 2 
with the Lifshitz transition again preserved in the pres- 
ence of weak interactions. It would be very interesting to 
extend the DMRG results presented here in order to ex- 
amine the possibility of a violation of logarithmic scaling 
at Lifshitz critical points with interactions present. We 
are currently investigating this issue. 

VII. CONCLUSION 

We have theoretically investigated the scaling of en- 
tanglement entropy across Lifshitz transitions in gapless 
fcrmions where the transition is driven by a change in 
the topology of the Fermi surface. In one-dimension, the 
prefactor of the leading contribution to the entanglement 
entropy in the long distance limit, which is known to ex- 
hibit a logarithmic correction, jumps as the number of 
Fermi points increases across Lifshitz transitions. The 
diverging length scale associated with the small Fermi 
pockets near the critical points shows up as a crossover 
scale in the entanglement entropy. In two-dimensions, it 



is shown that the leading and sub-leading coefficients of 
the entanglement entropy exhibit scaling behaviors with 
distinct exponents. Preliminary DMRG results for the 
one dimensional case indicate that similar behavior is 
present in weakly interacting systems. 
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